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ABSTRACT 

Using the STIS spectrograph on HST we have obtained a grid of [O III]AA4959,5007 and H/3 emission-line 
spectra at 0."05 x 0719" and 60 km s -1 (FWHM) resolution that covers much of the NLR of NGC 1068. We find 
emitting knots that have blueshifted radial velocities up to 3200 km s -1 relative to galaxy systemic, are 70-150 pc 
NE of the nucleus and up to 40 pc from the radio jet, emit several percent of the NLR line flux but no significant 
continuum, span a small fraction of the sky as seen from the nucleus, coincide with a region of enhanced IR 
coronal-line emission, show gradients in radial velocities of up to 2000 km s -1 in 7 pc, span velocity extents 
averaged over 0,"1 x 0,"2 regions of up to 1250 km s -1 , have ionization parameter U >0.1, and ionized masses 
~ 200 Mq/iL ',4 Ou .4 = 10 4 cm -3 ). The brightest parts of the blueshifted knots are often kinematically contiguous 
with more massive clouds nearer the jet that are moving with velocities of < 1300 km s -1 relative to galaxy 
systemic. However, some knots at 175 — 275 radii appear as bright points in a broken shell of radius ~ 0755 (40 
pc) that is expanding at up to 1500 km s -1 , implying a dynamical age of ~ 1.3 x 10 4 yrs. Between 275-475 from 
the nucleus, emission is redshifted relative to systemic, a pattern that we interpret as gas in the galaxy disk being 
pushed away from us by the NE radio lobe. We argue that the blueshifted knots are ablata from disintegrating 
molecular clouds that are being photoionized by the AGN, and are being accelerated radiatively by the AGN or 
mechanically by the radio jet. In their kinematic properties, the knots resemble the associated absorbers seen 
projected on the UV continua of some AGN. 

Subject headings: Galaxies: Active — Galaxies: Individual (NGC 1068) — Galaxies: Nuclei — Galaxies: Jets 
— Galaxies: Kinematics and Dynamics — Galaxies: Seyfert 


1. INTRODUCTION 

Emission-line profiles that extend over several thousand 
km s -1 are hallmarks of activity in galaxy nuclei of all luminosi¬ 
ties. Some AGN also show absorption lines at UV rest wave¬ 
lengths, originating from warm gas located outside the region 
that generates the nonstellar continuum. This gas is invariably 
blueshifted over a range of several thousand km s -1 from the 
galaxy systemic velocity, yet often shows discrete components 
that each span several hundred km s -1 (Hamann et al. 1997). 
What process accelerates the clouds? 

Spatially complete spectral maps of nearby active galaxies 
can separate clouds that are agitated mechanically by the AGN 
from those that are simply illuminated by AGN photons (e.g. 
undisturbed gas in an “ionization cone”). In this paper we 
report on observations made with the Space Telescope Imag¬ 
ing Spectrograph (STIS) and its medium-resolution (M) grat¬ 
ings to map the “Narrow-Line Region" (NLR) of the Seyfert 
galaxy NGC 1068. In this NLR, Alloin et al. (1983) sampled a 
cloud complex that spans ~ 200 pc, is centered ~ 210 km s -1 
blueward of galaxy systemic velocity, and has ~ 1670 km s -1 
FWHM (~ 2200 km s -1 according to the HST FOS spectra of 


Caganoff et al. (1991)) velocity extent. The complex of clouds 
discovered by Walker (1968), and mapped in detail by Cecil, 
Bland, & Tully (1990, CBT hereafter) and Arribas, Mediavilla, 
& Garcia-Lorenzo (1996) (see also Baldwin, Wilson, & Whit¬ 
tle (1987)), was shown from FOS spectra (Rraemer, Ruiz, & 
Crenshaw 1998) to be mostly photoionized by nuclear radia¬ 
tion, and to have bulk motion consistent with acceleration that 
is independent of radius (Crenshaw & Rraemer 2000b, CKb 
hereafter). A recent STIS spectrum (Rraemer & Crenshaw 
2000, KC1 hereafter) shows that emission-line centroids near 
the “continuum hotspot” (a feature with substantial scattered 
nuclear light located ~ 0717 [12 pc] N of the optically obscured 
nucleus) exhibit a trend towards greater blueshifted velocities 
as the ionization potential of the line species increases, a result 
that extends the earlier finding by Marconi et al. (1996) who 
found that coronal, near-IR lines are extended across a < 4" 
region. 

Section 2 describes the acquisition and reduction of the STIS 
spectra. We show in §3 that the fastest knots have been accel¬ 
erated to radial velocities blueshifted from galaxy systemic by 
>3200 km s -1 , a pattern not evident in the single STIS spec- 
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trum discussed by KC1 and CKb. We consider various mech¬ 
anisms capable of accelerating clouds and knots in §4, outline 
why the knots may be related to quasar “associated absorbers" 
in §5, briefly consider future observations to test this hypothe¬ 
sis in §6, and summarize our conclusions in §7. In this paper, 
we assume a distance to NGC 1068 of 14.4 Mpc, so \"= 70 pc. 

2. STIS LONG-SLIT SPECTRA 

2.1. Data Acquisition 

Deep spectra at high spatial resolution are necessary to map 
and characterize the various dynamical subsystems in this NLR. 
We have therefore used the STIS, splitting 14 HST orbits evenly 
between M-grating (R ~ 5,000) maps of the [O III] and H/3 
emission-line profiles, and less extensive L-grating (R ~ 1,000) 
maps of far- and near-UV line fluxes. We used the 0"2-wide slit 
(0."19 on sky) as a reasonable compromise between resolution 
and mapping efficiency. To maximize sensitivity, we binned 
pairs of pixels along the axis of wavelength dispersion dur¬ 
ing CCD readout of the M-grating spectra, yielding 60 km s -1 
FWHM resolution. 

We planned to obtain spectra both along and perpendicular 
to the radio jet. However, available guidestars restricted us to 
RA. 38°, close to the axis of the large-scale radio jet. After a 
standard STIS peakup on the nucleus, we performed a small- 
angle offset perpendicular to the slit for each map position. 
We obtained five long-slit spectra before a safe event delayed 
the next spacecraft visit sufficiently that the large changes in 
the orbit of HST prevented us from reusing target offsets from 
guidestars. This required us to spend an unanticipated half-orbit 
a year later to reacquire guidestars. The added overhead, to¬ 
gether with a slight positional shift between the two pointings, 
meant that we obtained a total of six spectra with parallel slits 
that cover ~2/3 of the area of the NLR, with spatial resolution 
0."05 along the slit (and jet axis) and 0"19 in the perpendic¬ 
ular direction, see Fig. 1. We supplemented these data with 
several G430L and G430M spectra from the HST archive (P.I. 
R. Antonucci) that were unbinned in wavelength and used the 
0."l-wide slit; the three archival M-grating exposures yielded 
a single spectrum along P.A. 10° that intersects the continuum 
hotspot, close to the P.A. 22° orientation used by a published 
STIS spectrum (KC1). Standard STIS processing with on-the- 
fly recalibration using optimal files produced wavelength and 
flux calibrated spectra whose parameters are listed in Table 1. 

2.2. Data Reduction 

We used the IRAF ‘cosmic’ task to remove charge spots and 
hot pixels. To derive line fluxes, we used the average red¬ 
dening derived by KC1 for the blue-wing of the line profiles, 
£g_y = 0.35 as an initial guess to be refined at different points 
along the slits during modeling of our more dust-sensitive UV 
spectra (Groves et al. 2001). This reddening is consistent with 
estimates from the [Fe II]A1.257 pm to 1.64 /on flux ratios 
across the NLR (Marconi et al. 1996). 

The spectra were examined for instrumental scattered light 
following a procedure similar to that of Nelson et al. (2000); 
none was found. After masking emission lines, we used the 
IRAF ‘continuum’ task with third-order Chebyshev polynomi¬ 
als to model and subtract continuum light on a row by row basis. 

From ground-based spectral maps (for example CBT, and 
Arribas, Mediavilla, & Garcia-Lorenzo 1996), line emission 
in this NLR is known to be structured into roughly a half- 
dozen major subsystems with different kinematics, excita¬ 


tion, and geometric distribution. The gas distribution appears 
to be even more complex in the HST Faint Object Camera 
(FOC) image (Fig. 1), where many knots are unresolved spa¬ 
tially. This implies that there will be a “slit effect” on de¬ 
rived radial velocities, the magnitude of which depends on 
where the knot falls within our 0" 19-wide slit. This effect 
is small compared to the line widths of features in our data, 
and amount to ±62 km s -1 (±1 resolution element) per max¬ 
imum ±0"095 displacement (STIS Instrument Handbook, at 
http://www.stsci.edu/instruments/stis). 

The lines of the [O III]AA4959,5007 doublet are often 
blended because of the large range of gas velocities. To clean 
line profiles for analysis, we therefore used two complementary 
approaches: 1) parameterize each line with up to 25 Gaussians 
and constrain the two lines to have the correct theoretical ratio 
where they overlap, 2) multiply the profile of the red wing of 
[O III] A5007 by the theoretical A4959/5007 flux ratio and then 
subtract the result from the observed red wing of the A4959 line. 
We synthesized the [O III]-only image from our spectral grid, 
and then aligned this image to the FOC image to better than 
0."1 accuracy (half of the slitwidth); we thus verified that the 
FOC F509N filter passes all velocity components in our spec¬ 
tra. Capetti, Macchetto, & Lattanzi (1997) have performed ab¬ 
solute astrometry and aligned the FOC image with VLA radio 
images, finding that radio feature S1 coincides with the warped 
maser disk, presumably the AGN, while S2 is a knot in the jet. 

3. EMPIRICAL RESULTS 

3.1. Kinematic Components in the [O III] and H/3 
Emission-Line Profiles 

3.1.1. Overview 

Our STIS spectra resolve the spatio-kinematic structure of 
this NLR, subdividing the major complexes into three major 
kinematic components that have different velocity substructure 
yet similar [O III] and H/3 line profiles: 

1. that associated with the rotating galaxy and its ioniza¬ 
tion cone (the narrow “spike” in the NE quadrant, see 
Alloin et al. 1983), 

2. redshifted filaments, mostly associated with the NE ra¬ 
dio lobe, which arise as the lobe pushes down on the 
disk, as originally suggested by Pecontal et al. (1997, 
PFBW hereafter), 

3. blueshifted components that, at fainter levels, extend 
to velocities of more than 3000 km s -1 from systemic 
(larger than the 2950 km s -1 separation of the [O III] 
doublet). 

We have associated kinematic components with individual 
clouds or cloud complexes located in the FOC image (Fig. 1). 
Slit numbers are indicated in Fig. 1 and also refer to panel 
numbers in the spectrograms of the seven slit positions Fig. 2. 
We focus on the various blueshifted components, summariz¬ 
ing their line profiles and the velocity dependences of their 
[O III]/H/3 flux ratios in Fig. 3. All velocities are quoted rel¬ 
ative to the systemic value of the galaxy (1148 km s -1 . Brinks 
et al. 1997). 
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3.1.2. Components HV0-HV7 

These features are isolated knots or diffuse edges of bright 
clouds in the FOC [O III] image, and are labelled 0-7 in Fig. 1. 
Their line profiles in the left columns of Fig. 3 show: 

1. A component whose flux centroid redshifts from sys¬ 
temic velocity at knot 0 to 600 km s -1 at knot 7. It has 
two FWHM values: in knots 0-3 it is 300-400 km s _I , 
while in knots 4-7 it is 900 km s -1 . The descriptions 
that follow omit this component. 

2. An extended blue wing out to -3200 km s _1 , whose be¬ 
havior in the different gas complexes we now detail. We 
refer to these as “high velocity" features HV0-HV7. 

Complex 0 is sampled by slits 1 & 2. In the FOC image, it 
is the diffuse, SE, boundary of the brighter, central NLR. Flux 
peaks near systemic velocity in the SW part of this complex and 
then extends to -2700 km s -1 as the distance from the nucleus 
increases by 072 (14 pc). 

Complex 1, centered 075 (35 pc) NE of complex 0, is a re¬ 
markable cluster of roughly a dozen unresolved knots in the 
FOC image. The knots subdivide into a comparable number of 
bright, blueshifted peaks in each of slits 1 to 3, and we can usu¬ 
ally match a knot to a velocity component because the knots are 
separated along the slit. The spectrum of the knot in slit 1 shows 
a diagonal feature, which is resolved into a series of spatially 
resolved steps in space-velocity that each span —400 km s -1 , 
indicating an increase in gas velocity with position along the 
slit from ~ -400 km s _1 to — -2500 km s -1 over 075 (35 pc). 
In slits 2 and 3, two increases with nuclear distance occur: from 
-250 to -1250 km s -1 then -1200 to -3200 km s -1 over the 0."4 
(28 pc) extent of the complex. Despite multiple discrete spatial 
components, the overall emission-line profile is smooth. 

Knot 2 appears only in slit 4, where it is distributed between 
-1800 and -300 km s -1 . It is elongated slightly N-S in the FOC 
image. 

Knot 3 is a spatially diffuse knot in slit 4, with gas between 
-100 and -2400 km s _l . The trend along slit 4, through knots 
2 and 3 and then farther NE, is that of an incomplete “Doppler 
ellipsoid" centered at -600 km s -1 and extending around that 
velocity by ±1500 km s -1 . In the FOC image (Fig. 1) it ap¬ 
pears as an elliptical ring of diameter 0755 centered 272 from 
S1 (50 pc) delineated by a half-dozen faint wisps (see pink bro¬ 
ken ellipse in Fig. 1). This apparently disintegrating structure 
therefore has a dynamical age of - 1.3 x 10 4 yrs. 

Complex 4 is included in slit 5, has 3 or 4 unresolved knots 
in the FOC image, and lies along the edge of complex G. Ve¬ 
locity features are evident in slit 5, and extend in two adjacent 
bands, from -200 to -2900 km s 1 and from -1200 (faint) to < 
-3000 km s -1 . 

Complex 5 forms the diffuse, NE boundary of clouds G and 
H. We may just resolve it in slit 6 into a linear gradient from 
-200 to -2300 km s -1 as nuclear distance increases. 


Knots 6 & 7 lie outside the area mapped by our six parallel 
spectra, but are intersected by the archival observation with slit 
along RA. 10°. These knots resemble others in the FOC image, 
with knot 6 extending to -1800 km s -1 and knot 7 extending to 
-1600 km s _1 . 

3.1.3. Clouds 

Our spectra in the right-hand columns in Fig. 3 include all 
the large clouds first delineated by Evans et al. (1991). Several 
clouds were spanned by CKa with their single L-grating slit, 
so the local continuum properties that those authors established 
are also mentioned in the following remarks. 

A, included in slit 3, has a small peak centered at 500 km s -1 , 
in addition to weaker emission that spans -1900 to +1000 
km s -1 . 

B, included in slit 4, extends from < -2200 to >+2500 
km s -1 (the velocity extent is uncertain because of the difficult 
continuum subtraction in our narrow spectral bandpass). The 
emission extends ~ 072 along the slit. The associated contin¬ 
uum is the “hotspot” discussed in detail by CKa that is probably 
scattered nuclear light. 

C, included in slit 5, extends from -1400 to +2000 km s -1 . 
CKa found that most of its continuum is scattered nuclear light. 
Bicknell et al. (1998) propose that this cloud is being agitated 
by knot C in the radio jet, which Roy et al. (2001) argue is a 
shock standing in the jet, based on its low proper motion in 
multi-epoch radio VLBI images. 

D is also included in slit 5. A faint unresolved component 
reaches -1800 km s -1 and attaches to brighter, spatially ex¬ 
tended emission that spans -500 to +1800 km s -1 . CKa found 
that half of its continuum is scattered nuclear light. 

E, included in slit 6, spans -1500 to +1300 km s -1 , but is 
brightest between -500 and +100 km s -1 . 

F, included in slit 3, has two separate plumes. That nearer 
to the nucleus extends from systemic velocity to -2700 km s -1 , 
that farther extends from -4-00 to +2000 km s -1 . 

G, included in slit 5, starts abruptly at -200 km s -1 and ex¬ 
tends to +1400 km s -1 . It emits strong [Fe VII] (PFBW). At 
lower spectral resolution, PFBW and KC1 both found a very 
large flux ratio [O III]/H/3 = 22 near peak emission. However, 
our high-resolution spectra isolates this component more effec¬ 
tively, reducing the summed flux ratio to ~ 17. 

H is sampled by slit 6. Its profile resembles that of cloud G, 
but is skewed to brighter emission at blueshifts between -500 
and -1300 km s~ 1 vs. -250 to -900 km s -1 for G. CKa found 
that half of its continuum is scattered nuclear light. 

J & K are covered by slits 4+5 and 6, respectively. They lie 
on straight lines drawn back to the nucleus from, respectively, 
clouds G and H. Axon et al. (1997) discovered that these fea¬ 
tures are double peaked, and we find that they are centered first 
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at +500 km s 1 then at -700 km s 1 as distance from the nucleus 
increases. 

I is the brightest line-emitting part of gas associated with 
the SW radio lobe. Line profiles are double-peaked, with the 
stronger peak redshifted (near +700 km s -1 ) and sharp in slits 
1-3, but more diffuse in slits 4-6. As the nuclear distance 
increases by 1", velocity centroids shift from +250 to +800 
km s -1 in slit 5 and from systemic to -500 km s -1 in slit 6. 

3.1.4. Clouds Associated with the NE Radio Lobe 

Our slits sample many filaments across the base of this ex¬ 
tended bow shock (Wilson & Ulvestad 1987) in the interval 2- 
4" from the nucleus. Some filaments beyond clouds G and H 
and a few fainter knots beyond the top of Fig. 1 are blueshifted, 
but beginning at 2 "6 radius most of what appears near the top of 
the FOC image as a coherent V-shaped feature is redshifted in 
slits 1-6, confirming the findings of PFBW. Between 2"9 and 
4" radius, it is striking that individual components are resolved 
along the slit yet all range in velocity from systemic to +1200 
km s -1 . There is a narrow feature near systemic velocity and 
only faint, diffuse emission to the blue. From the [S II] dou¬ 
blet in this region, PFBW derived gas densities near systemic 
velocity of ~ 400 cm -3 if T e = 10 4 K. 

3.2. Physical Properties 

We have used the line profiles displayed in Fig. 3 to constrain 
the ionized gas mass, kinetic energy, and momentum of each 
velocity component. A mask was used to isolate the component 
in the higher signal/noise ratio [O III] composite profile. This 
mask was then applied to the same velocity range in the dered- 
dened H i profile. The resulting flux has been converted into an 
ionized gas mass, assuming case-B recombination conditions 
(for which KC1 found evidence in the red-wing component NE 
of the nucleus) and a pure H plasma. This conversion depends 
on the average electron density. In Table 2 we have therefore 
scaled some quantities by density n e< 4 = 10 4 cm -3 , which Do- 
pita (1997) found was appropriate for the NLR. Undoubtedly, a 
large range of gas densities are present in each emitting feature. 

The high-velocity features HV0-HV7 emit a few percent of 
the total NLR H (3 flux. In Table 2 we also tabulate properties 
of clouds A-H that were identified by Evans et al. (1991), and 
introduce clouds J and K (cloud I is rather ill-defined, in the S W 
radio lobe). 

4. DISCUSSION 

While the clouds might have been as blended at HST reso¬ 
lution as they are from the ground, in §3 we showed that our 
deep spectra have resolved spatially the NLR. We have also 
resolved much of it kinematically, although not down to the 
thermal width (17 km s -1 in 10 4 K gas) in part because of inad¬ 
equate spectral resolution and in part because several compo¬ 
nents are still projected along each line of sight and may have 
turbulent substructure. 

To track the energy flow and to isolate possible shock-excited 
structures, we must now constrain the internal properties and 
space velocities of the NLR clouds. 

4.1. Photoionization of the NLR Clouds 

The NLR clouds are bright features embedded in the large- 
scale ionization cones and are distributed over several arc- 
seconds, so have been studied extensively from space and, in 


aggregate, from the ground. FOS spectroscopy lacked spatial 
resolution to separate clouds from their surroundings. Subse¬ 
quently, KC obtained a STIS spectrum along one cut through 
the NLR using low-resolution gratings that spans AA115 - 1027 
nm. Along the slit, they analyzed continuum fluxes (CKa), line 
fluxes (KC1), and gas kinematics (CKb). In KC1 they repro¬ 
duced the observed flux ratios with AGN photoionization mod¬ 
els that require considerable absorption along sight-lines from 
the nucleus that intersect the galaxy gas disk. They confirmed a 
correlation between ionization potential and line blueshift from 
galaxy systemic velocity that Marconi et al. (1996) had discov¬ 
ered in their ground-based IR spectra. 

KC1 sampled several clouds NE of the nucleus, so their 
photoionization models also address the new features that we 
have mapped. Their models combine matter- and ionization- 
bounded clouds, the former arising plausibly from low-density, 
photoevaporated envelopes of the latter (Binette, Wilson, & 
Storchi-Bergmann 1996). Matter-bounded clouds may filter 
AGN photons that impinge on the ionization-bounded cores, so 
internal dust would play a critical role in setting their emission¬ 
line flux ratios (Dopita et al. 2001; Groves et al. 2001). 

The models of KC1 required that the blue-wing clouds be 
optically thin above the Lyman limit, and need another pho¬ 
ton source between 1"4-1"8 to explain strong NUV lines and 
a larger He IIA4686/H/3 ratio along 202° RA. CKb suggested 
that clouds are decelerating here, forming a ~ 1000 km s -1 
shock that would brighten high-ionization lines as observed. 
They lacked the sensitivity and spatial coverage of Marconi 
et al. (1996) to coronal-line emission, so unlike those authors 
found no such emission beyond cloud B and the UV continuum 
hotspot. In CKa they established that the scatterer is tenuous 
and coincides with the deceleration region. They found reason¬ 
able fits to blue-wing line fluxes if 85-90% of the flux comes 
from tenuous gas with (logTY, Uh, log Nh) = (-1.45 ±0.15,3 x 
10 3 cm -3 ,21.18 ±0.18cm -2 ) and the rest from denser clouds 
with (log 77, n H , log Ah) = (-2.9,2 x 10 4 cm -3 ,20.7 ±0.1 cm -2 ), 
both components needing dust fractions of 10-25% (larger val¬ 
ues closer to nucleus) to suppress Lyo and C IV to their ob¬ 
served dereddened fluxes. 

KC1 found that the red-wing clouds have the higher dust 
fraction (25%) and are exposed to a lower ionizing flux (es¬ 
pecially above the He II Lyman limit), indicating absorption of 
AGN photons on the way to the red wing gas, and thence by the 
tenuous blue-wing component along our line of sight. Marconi 
et al. (1996) find IR coronal lines across the inner 4" diameter 
that emit out to 2000 km s -1 blueshift from systemic velocity, 
and which can arise either in gas photoionized by a hard UV 
continuum or in a hot (~ 10 6 K) collisionally ionized (perhaps 
shocked) plasma. Oliva et al. (2001) note that the [Fe II] emis¬ 
sion is unlikely to come from shocks because they find that the 
flux ratio [Fe II]/[P II]~ 1, 5% of the value seen in supernova 
remnants. 

X-ray results also support photoionization NE of the nucleus, 
beyond radius r,-. First, Young et al. (2001) show that fila¬ 
ments in the Chandra ACIS image correlate well with those in 
[O III] A5007 and Ha+[N II] images. Second, the Newton/XMM 
RGS spectrum reported by Paerels et al (2000) (also Behar et 
al. 2001) show that gas in this region is photoionized, based on 
the detection of narrow radiative recombination continua, line 
flux ratios within the helium-like triplets, and the weakness of 
Fe L-shell compared to K-shell emissions. 
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4.2. Deprojected Velocity Field 

Ground-based Fabry-Perot (CBT) and integral-field spectra 
(PFBW) have constrained the bulk motion and space distri¬ 
bution of the NLR clouds. CBT found an overall distribution 
of clouds in a thick bicone with maximum half opening angle 
~ 41° and whose axis is inclined ~ 85° to our line of sight; 
PFBW showed that gas coincident with the NE radio lobe be¬ 
yond 2." 5 radius is predominantly redshifted. CKb also posited 
a biconical flow to match the velocity centroids along their 
single STIS slit, requiring outflow velocities that reach 1300 
km s _1 at r, ~ l."7 NE of the nucleus and then decrease be¬ 
yond. Limited spatial coverage of a structure whose symmetry 
axis lies near the plane of the sky meant that their kinematical 
data are consistent with either an outflow that is radial from the 
nucleus or one that is roughly perpendicular to the jet axis. The 
archival M-grating spectrum (Fig. 4) shows that knots HV6 and 
7 are comparable to the largest blueshifts in the model geometry 
of CKb. 

We see abrupt jumps in velocity as the STIS slits cross the 
clouds and knots evident in the FOC narrow-band image. In 
Fig. 5, we compare our line profiles with the projections of sev¬ 
eral models of radial and cylindrical outflow. By increasing the 
half angle of the NE bicone from ~40 to 50° within the radius r, 
of the CKb model where outflow velocities begin to decrease, 
we better match the trend in the centroids of radial velocities 
of the blueshifted features in our spectra. The velocity limits 
of our two models are shown in the Figure, and correspond to: 
(solid line) increasing the maximum outflow to at least 2500 
km s -1 at r,, and (long dashed line) having the gas expand from 
the jet axis with all of its motion along our line of sight. We also 
plot the outflow model of CKb in short dashed lines. The gas 
becomes increasingly blueshifted (on average) out to radius 2", 
implying its acceleration on that scale. However, the observed 
gas does not fill all of the NE bicone, because redshifted emis¬ 
sion is observed at much lower velocities than predicted by all 
of the outflow models until the NE radio lobe is reached. Emis¬ 
sion then jumps to a range that extends from zero to > 1000 
km s' 1 redshift (relative to systemic velocity) and persists from 
3" to the middle of the NE radio lobe (which occurs at the top of 
the panels in Fig. 2); the emitting features are spatially compact. 
KC1 noted that the emission-line flux ratios of the red wing gas 
are consistent with irradiation by a heavily absorbed compo¬ 
nent. If these clouds are accelerated radiatively, they may not 
attain velocities as high as those that were irradiated directly. 

As suggested by Tecza, Thatte, & Maiolino (2001) based on 
infrared [Fe II] spectra and PFBW based on spatially decon¬ 
volved optical spectra, the red wing gas beyond 2 "5 radius is 
apparently being pushed into the dense galactic disk by the lat¬ 
eral expansion of the NE radio lobe. Our favored scenario is 
shown in Fig. 6. The jet is inclined to the disk. Eventually, at 
radius 2", it bursts out of the disk and inflates the radio lobe. 
The lobe expands laterally in all directions, but is only visible 
where it pushes into the disk and creates the zero to redshifted 
emission. The large blueshifted extensions on the high-velocity 
clouds suggest that the clouds are moving slower than the HV 
knots. The HV flow must have been accelerated relative to the 
clouds or reoriented closer to our line of sight by a combination 
of bouyancy, radiative acceleration, and lobe/cloud collisions 
(Taylor, Dyson, & Axon 1992). Depending upon their dynam¬ 
ics, the larger clouds may lag the knots because of larger drag 
forces, as KC1 have proposed for the red wing component. 


4.3. Acceleration of Clouds & Knots 
4.3.1. Jet Interactions with Clouds C & D 

Using the long-slit mode of the FOC, Axon et al. (1997) 
found that line profiles are double-peaked and have highest 
gaseous excitation across the radio jet near clouds G and what 
we call K. On this basis they argued that hot shocked gas in 
the jet cocoon expands clouds from the jet axis. Nearby, Krae- 
mer, Ruiz, & Crenshaw (1998) found evidence from emission¬ 
line ratios in an FOS spectrum for cosmic-ray heating of the 
emission-line gas. The Merlin and VLA images of Gallimore 
et al. (1996, G96 hereafter), when registered using the astrom¬ 
etry of Capetti, Macchetto, & Lattanzi (1997) (see red contours 
in Fig. 1), place NLR clouds B, C, D, and F adjacent to the ra¬ 
diojet; in particular, clouds B and C straddle a stationary radio 
knot (Roy et al. 2001) where the jet bends ~ 23°. Bicknell et 
al. (1998) have modeled this complex as a jet-cloud interaction. 
What constraints do our data impose? 

KC1 found that cloud C is photoionized. However, decom¬ 
posing its profile in Fig. 3 into Gaussians isolates a com¬ 
ponent centered ~ 300 km s -1 blueward of systemic and a 
wing component centered 300 km s -1 to the red, both with 
dispersions ~ 500 km s -1 . The blue component has flux ra¬ 
tio [O III] A5007/H/3 ~ 16 vs. ~ 8 —10 for the red component. 
While the considerable line widths are consistent with shocks, 
the flux ratios are consistent with photoionization by the AGN 
with the red wing gas more absorbed. Thus, deflection of the 
jet at cloud C is not clearly connected to the excitation of the 
cloud. Cloud E, which is farther from the jet than cloud C, 
shows a similar blue feature but not one in its red wing. 

The NE component of the nuclear radio triple is 0."75 NE 
of knot SI (the nucleus). G96 note that this jet component, like 
the radio knot near cloud D, has a flat synchrotron spectrum and 
extended emission that is misaligned with the jet. However, the 
jet does not bend at the NE component nor is there maser emis¬ 
sion, arguing against a shock front. It is therefore interesting 
that the emission-line profiles of cloud D are double peaked. 
The blue peak has flux ratio [O III]/H/3 ~ 15, but now the red 
peak ratio decreases monotonically through the profile from 17 
to < 3 over a 1000 km s -1 range. This decline implies that gas 
density increases with the velocity deviation from galaxy sys¬ 
temic, not what would be expected if most of the cloud mass is 
at velocities close to rest in the galaxy. G96 note that a cloud 
interacting with the jet must be sufficiently massive to avoid 
acceleration, implying 

M ci ouiB > 40M o ( 7 ^-) 

f _ Yl _ ( tlobe ^ 

^ 15,000 kms~ l ) V3.5x 10 5 yrJ ' 

Here jet density and velocity are normalized to the values ap¬ 
propriate for the lighter jet model of Wilson & Ulvestad (1987) 
for the NE radio lobe, t/ 0 b e is the age of the NE radio lobe from 
spectral aging as estimated by G96, and the jet radius rj was 
measured by Gallimore et al. (1994). If the average NLR cloud 
density is ~150 cm -3 , the entries in Table 2 show that clouds are 
massive enough to remain adjacent to the deflecting jet. This 
value is consistent with the mean densities of Galactic molecu¬ 
lar clouds (Solomon et al. 1987). 

Elsewhere, radio and line emission are uncorrelated. G96 
note that this is expected if the jet has swept all gas from the 
region. The high-velocity knots might then be ablated from 
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molecular clouds in the galaxy disk that have rotated into the 
jet. The clouds accelerated would all be blueshifted because 
they would move predominantly away from and above the 
galaxy disk. Similar clumps below the jet would sink into the 
denser disk where they would be subject to stronger drag forces 
and hence decelerated more effectively (consistent with the gas 
distribution in the spectral maps of PFBW). 

4.3.2. Shocks Near the X-ray Subpeak? 

Three characteristics of the XMM/Newton X-ray spectrum of 
NGC 1068 within 20" of the nucleus indicate that the gas is 
predominantly photoionized: (i) narrow radiative recombina¬ 
tion continua; (ii) in the He-like triplets the forbidden lines are 
stronger than the resonance lines; and (iii) Fe L emission is 
weak compared to the K shell emission of the abundant ele¬ 
ments (Kinkhabwala et al. 2001; Behar et al. 2001). However, 
subtraction of the predictions of a photoionization model from 
the spectrum shows residuals that may indicate a collisionally 
excited component or other effects (Paerels et al 2000). 

Moreover, near the V-shaped line-emitting feature whose 
apex is 3" NE of the nucleus (at top of Fig. 1), the 4-26 A 
spectrum from the Chandra HETGS (Ogle et al. 2001) shows 
that resonance lines are brighter than forbidden lines, imply¬ 
ing either a collisionally ionized component or strong resonant 
scattering. Ogle et al. (2001) show that the collisional compo¬ 
nent lies at the X-ray subpeak reported by Young et al. (2001). 
The X-ray subpeak may represent gas in the galactic disk on the 
far side of the NE lobe and being compressed by it (cf. PFBW). 
As discussed in §4.1, photoionization models of the UV and 
optical STIS spectra of KC1 also support this additional power 
source near the subpeak. 

4.3.3. Shocks Near Clouds G-K? 

Two lines of evidence argue for shocks near clouds G-K: 

1. Clouds G and K are close to the “NE radio subpeak", 
located 2" 1 NE of the obscured nucleus (Wilson & Ul- 
vestad 1987), suggesting dissipation of kinetic energy 
and possibly cloud acceleration by the radio jet. 

2. Nearby we find the Doppler ellipsoid described in 
§3.1.2, whose expansion velocity of 1500 km s -1 will 
generate a post-shock temperature of ~ 3 x 10 7 K if the 
ambient medium is relatively stationary, fully ionized, 
and has 10% He abundance. 

To estimate the mechanical luminosity of the Doppler ellip¬ 
soid, we assume that energy is injected at constant rate, that the 
shell shock is strong so that we can ignore the thermal pressure 
of the ambient gas, and that the bubble is within the scale height 
of the galaxy disk so that the ambient ISM density is constant. 
Then in the energy conserving phase of the bubble interior, the 
bubble radius and expansion velocity are related: R = 1 .61 vt. 
By observation, t = 0.012 Myr. Mechanical luminosity L fol¬ 
lows from R(t) = 3\{L^,( s /nof , ■ 2 t®■ (, P c = 100 P c with { 6 = 0.012 
(where has units 10 36 ergs s -1 and no is the density of the 
ambient ISM in cm -3 ), so L = 5.2 x 10 7 «o Lq. Overall, the fil¬ 
aments that delineate the bubble and maybe its associated ISM 
shock emit <2% of the NLR [O III] and H (3 fluxes. 

4.3.4. Radiative Acceleration of Dusty Clouds 


Dust absorption dominates opacity in a photoionized plasma 
like this NLR when ionization parameter U > a(T e )/cn, where 
c is lightspeed, k ~ 10 -21 cm 2 atom -1 is the opacity assuming 
the standard Galactic reddening curve, and a(T e ) ~ 2 x 10 -13 
is the recombination coefficient. The dust will compete suc¬ 
cessfully for ionizing photons in the highly ionized zone, sup¬ 
pressing line flux. A strong pressure gradient will also be set 
up so that the plasma density will be much enhanced by the 
time lower ionization species emit. In consequence, a lower 
U ~ 0.007 will provide an apparent model fit to the spectrum, 
close to the value preferred by KC for their denser compo¬ 
nent. This ionization parameter is closely related to the value 
of U at which radiation pressure starts to dominate either the 
pressure gradient, or the dynamical acceleration of the ionized 
plasma. Thus, the flow of dusty clouds falling toward the nu¬ 
cleus is intrinsically unstable above the dust sublimation point, 
with clouds being driven into outflow in the direction of slowest 
accretion (i.e. within the ionization cones). Dopita et al. (2001) 
detail this scenario. In hydrodynamical simulations of cloud ac¬ 
celeration (Schiano, Christiansen, & Knerr 1995), tenuous gas 
streams ablate from surfaces of massive clouds, which may be 
the source of the high-velocity features we observe. 

In this NLR, the flux in the ionizing continuum needed to 
maintain the observed line luminosity (~ 3 x 10 44 ergs s *, 
KC1) is comparable to the FIR luminosity 1.5x 10 11 L 0 . Dusty 
clouds are accelerated radiatively by ~ 10 -5 cm s -2 , so can 
reach their observed mean velocities of ~ 500 km s -1 over 35 
pc. In Fig. 3, the cloud line profiles indeed show less velocity 
structure than do the knots. 

This scenario ignores the role of the radio jet in agitating the 
NLR clouds (Bicknell et al. 1998). But, as discussed earlier, its 
influence may be localized to line broadening and ablation in 
specific jet/cloud and lobe/cloud interactions, for example be¬ 
tween B and C (Gallimore et al. 1996) and G and H (Axon et 
al. 1997). While the region of collisional line emission coin¬ 
cides with the base of the NE radio lobe, numerical simulations 
of a breakout lobe do not show (D. Balsara, private commu¬ 
nication) a strong backflow vortex at its base that might form 
high-velocity shocks. In any event, many of the high-velocity 
knots have similar line profiles, and are found > 40 pc from the 
brighter parts of the jet. 

5. ARE THE HIGH-VELOCITY KNOTS “ASSOCIATED 
ABSORBERS” SEEN IN EMISSION? 

If the HV emission knots in NGC 1068 were instead viewed 
against the nuclear continuum, they would resemble in their 
kinematics the Associated Absorption Line (AAL) systems 
seen in a few percent of quasars (Hamann et al. 1997, for ex¬ 
ample). KC1 argue from their spectral models that the putative 
UV absorber in the NE quadrant must block less than half of 
the nuclear light. If the blueshifted knots have the same extent 
along the line of sight as we see on the sky (assumed here to be 
just resolved), they would cover at most 2% of the sky at the nu¬ 
cleus, consistent with values derived for some AAL absorbers. 

Such absorbers are usually thought to lie much closer to the 
source of hard photons than the 35-160 pc radial distance of the 
clumps in NGC 1068. However, recent photoionization mod¬ 
els of high-resolution quasar spectra are establishing that AAL 
absorbers lie at diverse distances from the AGN, including a 
range comparable to those we see. For example, de Kool et al. 
(2001) established that some absorbers in a BALQSO must be 
~700 pc from that nucleus and have kinematic structure sim- 
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ilar to the HV components in NGC 1068. High-velocity fea¬ 
tures are also visible on similar scales in STIS spectra of NGC 
4151 (Hutchings et al. 1999), although at more modest veloci¬ 
ties (1400 km s -1 ) than we see, and are not clearly connected to 
massive clouds as in NGC 1068. 

6. FUTURE WORK TO CONSTRAIN FILAMENT PROPERTIES 

The brighter blueshifted knots in NGC 1068 could be 
mapped in higher ionization lines with UV spectra to estab¬ 
lish their optical depths and ionization structure. But this study 
would be arduous with STIS because of their relative faint¬ 
ness and large reddening. Unfortunately, the complex back¬ 
ground emission makes the knots unsuitable targets for the 
high-efficiency Cosmic Origins Spectrograph (COS) to replace 
ST IS on I 1ST. 

Now that specific targets are known up to 2" from the nu¬ 
cleus, practical spectroscopy is possible with optical and near- 
IR lines using integral-field spectrometers that are appearing 
on 8-10m telescopes. An interesting constraint follows if such 
spectra find evidence of the H I/II recombination front present 
in optically thick clouds (Shields, Ferland, & Peterson 1995; 
Viegas-Aldovandi & Contini 1989, for example). Constraining 
the column will better constrain photoionization models to pre¬ 
dict the full mass — hence KE, momentum, and ultimately the 
origin — of each knot and cloud. 

7. SUMMARY 

Our STIS spectral maps have resolved the spatio-kinematic 
structure of the NLR in NGC 1068, and reveal that many com¬ 
pact knots in the FOC image have large, exclusively blue- 


shifted velocities, ranging up to ~ 3200 km s -1 from systemic. 
If these knots are optically thick in the UV continuum, they 
are good candidates for “associated absorber" clouds present in 
other AGN. Lying 70-150 pc from the nucleus, several form a 
broken spatio-kinematic ring of diameter 35 pc that is expand¬ 
ing at up to 1500 km s -1 . The ring is adjacent to the more 
massive but slower moving NLR clouds G-K, so may originate 
from ablata streams from disintegrating molecular clouds that 
are being photoionized and accelerated radiatively by the AGN 
or mechanically by a knot in the radio jet. 

By resolving this NLR, we have demonstrated that ground- 
based integral-field spectrometers will be able to constrain the 
properties of the emitting gas in NGC 1068 by using classical 
plasma diagnostics in the optical and near-IR. More fundamen¬ 
tally, we may have a convenient view onto gas associated with 
an important component of AGN that is hard to study elsewhere 
because it is usually seen in absorption. 

This work is based on observations made with the NASA 
Hubble Space Telescope, obtained at the Space Telescope Sci¬ 
ence Institute, which is operated by the Association of Univer¬ 
sities for Research in Astronomy, Inc., under NASA contract 
NAS5-26555. These observations are associated with proposal 
ID GO-7353. We thank Alessandro Capetti for the HST FOC 
mosaic of NGC 1068, and NASA for support through grants 
NAG 81027 and GO-7353. GC thanks the Research School 
of Astronomy and Astrophysics in the Institute for Advanced 
Studies of the Australian National University, and its Director 
Prof. J. Mould, for financial support and a stimulating work en¬ 
vironment to begin this paper. 
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Fig. 1.— FOC image of [O III]AA4959,5007 line emission across the circumnuclear region of NGC 1068 (Macchetto et al. 1994). The greyscale is proportional 
to the log of the intensity. The vertical direction is in P.A. 38°, with NE at top. denote the high-velocity regions (HV0-HV7) described in §3.1.2; they are also drawn 
in Fig. 2 to define the relationship between this image and the locations of the slits. The numbers 0-7 of these clouds are given outside the figure at left and right. 
The vertical scale (white line and vertical numbers) originates at the (obscured) nucleus (which coincides with radio source SI, denoted by +) and is in arcseconds. 
The numbers at top and the white 7 in the figure refer to the slits, the boundaries of which are shown as vertical and diagonal lines. The spectra from these 7 slit 
positions are shown in Fig. 2. The location of the Doppler ellipsoid is dashed in magenta. The V-shaped feature discussed in §3.1.4 and §4.3.2 is at top, between 3 
and A" radius. 
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Fig. 2.— Continuum-subtracted, log-scaled intensity, spatial-spectral images of [O III]A5007 emission-line profiles, with position along the slit vertical and 
wavelength horizontal. The A4959 profile has been subtracted. Slit P.A. is 38°, except for panel (7) where it is P.A. 10° (marked by the parallel lines in panel 8, 
which repeats the FOC image shown in Fig. 1). The archival spectrum in panel (7) has been binned spatially to resemble spectra in panels (1-6). The position of the 
slit w.r.t. the circles and ellipses denoting the high-velocity (HV) complexes (see Fig. 1 and Table 2) is shown in panels (1-6). The vertical scale is in arc-seconds 
from the nucleus (radio source SI) and the horizontal scale is in km s -1 relative to the galaxy systemic value. Note the very broad (>5,000 km s -1 FWZI) [O III] 
profiles at knot B in slit 4. The Doppler ellipsoid (§4.3.3), which is dashed in magenta in Fig. 1, has its kinematic manifestation marked by a black dashed oval in 
panels 2-6. A high-fidelity, color version of this figure is at http://www.thececils.org/science/nl068/spectra.jpg 
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Fig. 3.— Shown are the dereddened line profiles of the high-velocity knots (HV0-HV7: left diagram) and more massive clouds (A-K: right diagram). In each set 
of three spectra, the top profile (a) is [O III]A5007 after contamination by A4959 has been removed, the middle (b) is H (3, and the bottom (c) shows the range of the 
ratio [O III]A5007/H/3 at each velocity within =bl cr statistical plus estimated systematic uncertainties from the continuum subtraction. Velocities are relative to the 
galaxy systemic value. The spectral resolution is 60 km s -1 FWHM for the [O III] and H (3 profiles; the line ratio plots have been smoothed with a boxcar of width 
160 km s -1 . The vertical scale is in units of 10 -18 ergs s -1 cm -2 pixel -1 , where the pixel bin width is 33 km s -1 . 
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Radial Velocity (km/s) 

Fig. 4.— Velocities of the flux centroids of the two main emission-line components in each continuum-subtracted line profile as a function of distance from the 
nucleus, obtained from an archival M-grating spectrum (slit P.A. = 10°). NE is at top and SW is at bottom. The vertical dotted line shows the systemic velocity 
of the galaxy. The velocity field is similar to that shown in Fig. 2 of CKb that was extracted along a slightly different P.A. The high-velocity, blueshifted motions 
between lf'5-275 radius NE of the nucleus result from knots 6 and 7. 


Table 1 

STIS Long-Slit Observations of NGC 1068 


Data Set 

Date 

Grating 

Dispersion 
(A pixeD 1 ) 

Spectral 

Range 

Exposure 

(s) 

Slit P.A. 

(deg.) 

O56502010 

1999 Oct 2 

G430M 

0.6 

4818-5104 

2585 

38 

O56502020 

1999 Oct 2 

G430M 

0.6 

4818-5104 

2775 

38 

O56502030 

1999 Oct 2 

G430M 

0.6 

4818-5104 

2775 

38 

O56502040 

1999 Oct 2 

G430M 

0.6 

4818-5104 

2775 

38 

O56502050 

1999 Oct 2 

G430M 

0.6 

4818-5104 

2772 

38 

O56503010 

2000 Sept 22 

G430M 

0.6 

4818-5104 

2294 

38 

O56503020 

2000 Sept 22 

G430M 

0.6 

4818-5104 

2853 

38 

O56503030 

2000 Sept 22 

G230L 

1.6 

1568-3184 

2376 

38 

O56503040 

2000 Sept 22 

G230L 

1.6 

1568-3184 

2376 

38 

O56503050 

2000 Sept 22 

G230L 

1.6 

1568-3184 

2376 

38 

O56503060 

2000 Sept 22 

G230L 

1.6 

1568-3184 

2376 

38 

O56501010 

1999 Sept 22 

G140L 

0.6 

1140-1730 

1425 

38 

O56501020 

1999 Sept 22 

G140L 

0.6 

1140-1730 

1425 

38 

O56501030 

1999 Sept 22 

G140L 

0.6 

1140-1730 

1425 

38 

o51j01050 a 

2000 Jan 14 

G430M 

0.3 

4818-5104 

800 

10 

o51j01060 a 

2000 Jan 14 

G430M 

0.3 

4818-5104 

846 

10 

o51j01070 a 

2000 Jan 14 

G430M 

0.3 

4818-5104 

960 

10 

o51j01080 a 

2000 Jan 14 

G430L 

0.6 

2900-5700 

720 

10 

o51j01090 a 

2000 Jan 14 

G430L 

0.6 

2900-5700 

720 

10 

o51j010a0 a 

2000 Jan 14 

G430L 

0.6 

2900-5700 

821 

10 


Archival exposure. 
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Fig. 5.— Same display format and region as for panels 1-6 in Fig. 2, but now with contours to show the space-velocity projections of the outer cones in the three 
bi-conical outflow models described in the text: (solid lines) increasing the maximum outflow velocity from the nucleus to 2500 km s -1 at r*, (long dashed lines) 
having the gas expand cylindrically from the jet axis at close to the observed velocity, (short dashed lines) outflow model of CKb. At best these models match only 
the blue wing of the data within 1"5 NE of the nucleus. The red wing is not fit at all. 
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Fig. 6.— Kinematic model of the NLR refined from that proposed by Tecza, Thatte, & Maiolino (2001). To represent patterns seen in our spectral maps, we have 
(i) added radial expansion from the nucleus, the velocity of that outflow appearing to increase with its angle of travel above the disk; (ii) removed blueshifted gas 
associated with the expanding radio lobe because our data show only redshifted gas that is being pushed into the galaxy disk; (iii) added the gray arrow to represent 
possible ablata from the molecular ring (which may be a highly warped outer extension of the putative torus obscuring our view of the AGN). 


Table 2 

Properties of the NLR Clouds in NGC 1068 


Cloud 

Name 

Area a 

(arcsec 2 ) 

L[0 III] A5007 b 
(10 40 ergs s -1 ) 

Ionized H Mass 
(10 2 M Q /n e , 4 ) 

[O III] A5007 d 
/H/3 flux 

log KE “ 
(ergs/n c , 4 ) 

log momentum e 
(dynes s/n^) 

A 

0.03 

3.1 

16.6 

12.3 

51.50 

43.87 

B 

0.04 

10.8 

63.0 

11.5 

52.33 

44.61 

C 

0.06 

10.9 

55.3 

13.2 

51.98 

44.38 

D 

0.04 

21.1 

100.8 

13.9 

52.18 

44.64 

E 

0.09 

20.7 

96.9 

14.2 

52.16 

44.58 

F 

0.06 

24.5 

116.2 

13.9 

52.58 

44.88 

G 

0.06 

8.1 

32.5 

16.7 

52.02 

44.37 

H 

0.06 

3.7 

16.6 

14.7 

51.75 

44.11 

J 

0.06 

3.3 

17.6 

12.7 

51.53 

43.97 

K 

0.06 

3.5 

18.8 

12.4 

51.65 

44.07 

HV0 


3.6 

20.6 

10.9 

51.79 

44.09 

HV1 


5.0 

26.4 

11.8 

52.17 

44.35 

HV2 


1.4 

8.1 

11.3 

51.31 

43.65 

HV3 


0.6 

5.5 

7.0 

51.75 

43.84 

HV4 


0.6 

4.2 

9.8 

51.48 

43.61 

HV5 


0.8 

5.0 

11.2 

51.33 

43.60 

HV6 


1.1 

6.2 

12.2 

51.49 

43.74 

HV7 


2.1 

13.9 

10.3 

51.60 

43.96 


“Area of each distinct kinematic component on the sky, typically uncertain by ±0.015 arcsec 2 . 
b Line fluxes at velocities below galaxy systemic were dereddened by E(B-V)=0.35, above systemic by 0.22. 
“Assumes case-B recombination conditions and a pure H plasma. 
d Velocity averaged flux ratio, typically uncertain by ±0.3. 

“Minimum values, assuming that all velocity is along the line of sight. Density n f ,4 is in units of 10 4 cm 4 . 







